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•niis report presents the technical basis for the design of 

the Hydrogen-Air-Vitiated Heater. The heater liner is subjected 

to a maximum thermal environment at specified conditon z' 

(Figure I-l), where the combustion gas temperature, pressure and 

flow rate are 5000 F, 750 psia, and 11.0 Ib/sec, respectively, 

2 

euid results in a heat flux of the order of 275 BTU/sec-ft . Cool- 
ing and stress analyses indicate that water is the logical 
choice for cooling of the combustor liner. A cylindrical shell 
of zirconium copper was selected as the combustor liner. This 

material, a high copper alloy, was chosen primarily because of 

2 

its high thermal conductivity (k = 200 BTU/hr-ft ) as well as good 
yield strength (35,000 psi) in a forged condition. Additionally 
by using a water cooled liner there is built into the design, 
potential for future extensions to more severe thermal environ- 
ments beyond the present specification envelope. 

A mixing analysis was undertaken to establish a good com- 
bination of combustor length and injector configuration. The 
analysis, using a conservative analytical approach, indicates a 
combustor length of the order of 5 ft combined with discrete 
fuel and oxidizer injection at an approximate 2-1/2 inch radial 
combustor position, and results in uniform combustion products 
at the heater exit for all specified envelope conditions. 

Equilibrium composition and performance curves were pre- 
pared to permit rapid determination of air, oxygen and hydrogen 
gas flow requirements relative to total gas flow, as well as total 
flow relative to heater pressure, temperature, and nozzle throat. 
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The equilibrium composition of the gaseous system was computed 
utilizing existing computer programs described in detail in GASL 
Technical Report No. 676 (Ref. II-l) • This repxjrt presents the 
results of a recent analytical investigation of the effects of 
vitiated air contamination on combustion and hypersonic airbreath- 
ing engine ground tests. Shifting equilibrium calculations for 
hydrogen vitiated facility combustors and nozzles were performed 
by minimizing the Gibbs Free Energy of the entire system, subject 
to the constraint of element conservation, and computing the 
equilibrium composition of the reacting gas. In addition, finite 
rate chemistry calculations were made employing reaction mechanisms 
and high speed computational techniques. The results indicate that 
the assumption of chemical equilibrium is adequate for facility 
nozzle flow determination over the range of tunnel conditions of 

interest. 

The equations used to compute the combustion gas properties 
are based on the following assumptions: 

1. The process is adiabatic and one-dimensional; 

2. The gas is in chemical equilibrium; 

3. The nozzle expansion is isentropic. 

The initial species concentrations of hydrogen, oxygen, and 
nitrogen were computed from the following reactions: 


Y y Y Y *5 

-JL H + — ^ O + N - H-0 + 0.232 0 + ~ N 

2 2 32 2 28 2 18 2 2 28 2 


1 Ibm Air - 0.232 Ibm 0, + 0.768 Ibm N 
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where the product species have been evaluated at conditions 
corresponding to flow in the test section. Thus, there will be 
23.2 percent oxygen by weight in the combustion products in the 

wind tunnel test section. 

The flow is assumed to be steady and one-dimensional, and 

ideal gas is assumed throughout. Under these conditions, the 

relevant, ecjuations arcs 

Equation of States 

- PM 
P “ RT 


Momentum Conservation: 


pu 


dx dx 


Energy Conservation: 

dh. dX^ 

? ^ ”i ST 


_ du 

3 ^ ^ 


= 0 


Continuity: 

w = puA 
or 

d£ + ^ ^ = 0 

p U A 

The results of these calculations are presented in 
Figures II-l through II-5 and permit the rapid determination 
of air, oxygen, and hydrogen propellant requirements relative 
to total propellant flow, as well as total propellant flow 
relative to heater pressure, temperature, and nozzle throat area 
Curves of the reactant species concentrations for hydrogen 
vitiated air are shown in Figure II-l. The reactants were in- 
jected into the burner at an initial temperature of 100 F. The 


6 


mmnii imuD tMwa iAto«ATe>*n ik 



results indicate the compositions required to increase the 
temperature of "clean" air from ambient to stagnation temperatures 
in excess of 5000°F for stagnation pressures between 300 and 
1200 psia. At the lower pressures more energy is absorbed into 
the increased dissociation, resulting in lower stagnation 
temperatures. The maximum temperature attainable, with oxygen 
replenishment, is determined by the case of zero air injection. 

The product species resulting from this vitiation, exp.xnded 
to conditions corresponding to flow in the facility test section 
where the primary species are frozen in composition, are shown 
in Figure II-2. These results arc for a gas in chemical equil- 
ibrium, and are plotted as a function of tunnel stagnation 
temperature. The molecular weight variation of hydrogen vitiated 
air is also shown. It is observed that increased vitiation 
(higher fuel injection) acts to decrease the molecular weight of 
hydrogen vitiated air due to the increased concentration of low 
molecular weight hydrogen. 

Figure II-3 presents a plot of the vitiated air mass flow 
parameter (.Ij/P^A*) as a function of tunnel stagnation temperature. 
This curve permits rapid determination of the total mass flow 
rate of propellant. Figures H-4 and II-5 present curves of the 
mass flow parameter specialized for a particular nozzle having 
a throat area equal to 2.38 in^ and operating at stagnation 
pressures of P^ = 1103, 804, 550, 294 psia. 
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III. MIXING ANALYSIS 


K. CENTRAL AND RADIAL CONFIGURATIONS 

In performing a mixing analysis for optimizing the heater design 
there are initially two chief governing criteria. The first is 
to ensure that the ignitor system will cause ignition of the 
hydrogen-oxygen-air flows in a short distance, and the second is 
to achieve a uniform hot flow in as short a heater length as 
possible. Since the heater will be started up in a gradual, 
transient manner there is no problem envolved in igniting all of 
the fuel. Hence, we are solely concerned with minimizing the 
length required to obtain a desired degree of uniformity in the 

flow field. 

The basic mixing analysis employing the boundary layer 
parabolic equations is described in detail in Refs.III-1,2.3, and 
has been applied to a wide variety of problems at GASL. Briefly, 
the conservation equations for a steady-state axisymmetric flow 
system with equilibrium burning may be written: 

Continuity: 

^(PUYl == 0 

&x 3y 


Momentiim: 


du du 


dE 1 

" dx y ^y 


IS 4 . ^ 


V> 2S .i . 1 


1 a t r 1 _ M 4 . f 1 - i = + Z 

y pr dy ^ pr ^ ay i.i 




— ) 3 } 

Pr^ by •* ^ 


Diffusion! 


daj b(yj 
(a) elements: pu _ + pv ^ 


1 a_ . jai . 

y dy ' Sc by 


(b) species; 


Pu 


_ + pv — = w + y ay ( sc ay 


The von Mises 


stream function transformations 


<l>y = p«y 


^ = p'^y 


are used to transform the 90 veming equations to a streamline 
coordinate system where the continuity equation is now explicit 
in the definition of the stream function: 

Momentum; 

( a— ) 

dx pu dx 


Energy: 


r SH 


M = i ^ r ^ + (Pr 
Sx fj) pt a<j) 


- 1) 


»( r > 

a^ 


^ i Iffl n 1 

+ r h (Le - 1) ^ J J 

i=l ^ 


Diffusion; 


ax 


1 ) 
^ ^ Sc 


where 


and 


a = 


_ upuy 


^8 «» ^ i i 

H- 

and the chemical system employed^ 
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Using a finite difference technique, as described in Ref JII-1. 
it is possible to obtain a numerical solution of the above set of 
equations providing that adequate models for the turbulent eddy 
viscosity, 4 and the chemical burning process are provided. 

Semi-empirical turbulent mixing models have been developed 
for ducted flows as described in Refs. 111-3,4.. Since we are 
concerned with the mixing problem, and the residence time in 
the low speed heater of a typical mass element is far in excess 
of the time required for hydrogen-air mixtures to burn completely 
and reach chemical equilibrium, a simple "complete combustion" 
chemistry model was employed. The difference between the complete 
combustion chemistry model and equilibrium chemistry is that the 
possible dissociation of and H^O into H,0, and OH at high 

temperatures is ignored. Thus, the final calculated temperature 
level of the flow is somewhat higher than it is in physical reality. 
However, this does not significantly effect the mixing process as 
will be shown by a comparison of an equilibrium chemistry calculation 

with a coinplGt© conibustion calculation* 

The calculations were performed at the extremes of the desired 
performance conditions of the heater. The initial conditions of 
each calculation are presented in Table III-l. 

The philosophy adopted in performing the analytical calculations 
was to try to make each step as conservative as possible. Thus 
the first heater fuel injection configuration to be considered 
was to have concentric rings of hydrogen, oxygen, and air about 
the axis. Prom the viewpoint of mixing, this is the worse possible 
configuration since the fuel has a minimal surface area in contact 
with the oxidizer. 

Constant area calculations were begun for the above 
configuration and in all four cases local flow reversal (velocities 
in the upstream direction) were encountered. This phenomena is 
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partly due to the very low initial velocities in the air stream. 
One-dimensional calculations showed that there was no possibility 
of thermal choking occurring in the heater. Hence, the initial 
flow reversal was a local eddy phenomena. 

A boundary layer type analysis cannot be used to analyze 
a flow field where local flow reversal occurs. However, by 
specifying that the flow be considered constant pressure instead of 
constant area, this difficulty was bypassed. Since the flow in the 
heater is quite subsonic (bulk Mach numbers varied from .01 to.l) 
the variation in static pressure between the constant area 
calculations and constant pressure calculations was less than 1% in 
three of the cases, and a maximum of 5% in case E . 

It should also be noted that having local flow reversal enhances 
the mixing process, and hence makes the calculations even more 
conservative. As is shown in Figs.III-1 through 10, satisfactory 
burning and mixing is achieved in a five foot length for high temperatur 
cases D and E. However for the low final temperature cases F and H, 
as much as twelve feet is needed to get a reasonable degree of 
uniformity. 

Upon obtaining this result, a more realistic configuration was 
considered. Here, the ring of hydrogen injectors was located at 
a radial distance of three inches from the axis, with the ring of 
oxygen injectors immediately above it. Complete combustion 
calculations were performed for cases D and F , since they represented 
the best and worse mixing performance, respectively, for the central 
configuration. 
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As is shown in Fig.III-H through 15, satisfactory uniformity 
is obtained for case F using the radial configuration, and the 
mixing performance of high temperature case D is not adversely 
effected by the change from the central to the radial configuration. 

Thus, the above calculations indicate that it is possible, 
using the radial configuration of hydrogen and oxygen tubes, to 
design the heater to provide a satisfactory degree of mixing 
within a five foot length, over the entire range of desired 
operating conditions • 

The final calculation to be performed involved repeating tn 
radial configuration calculation for case F, using equilibrium 
chemistry instead of the complete combustion chemistry model. When 
the results shown in Figs.IH-16 and 17 are compared with Figures 
ni-12 &14 it is clear that the use of the complete combustion 
chemistry model has been justified, since the more physically 
realistic equilibrium chemistry model indicates slightly better 
mixing. Thus, the complete combustion model, as desired, is 
conservative, and may be safely used for the design of the heater. 

One final point to be considered is what the effect would be 
if it were necessary to increase the velocities (for purposes 
of improving fuel injection system control and stability) at 
which the hydrogen and oxygen flows enter the heater. 

ZalOcay, in Ref. Ill -5 found experimentally that the potential 
core length of a jet (defined as the length required for a one 
percent change in chemical composition on the centerline of the 
jet) was determined by 

*o ■ ’'l 'jet '^(pu) outer ‘ 
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where is constant. 


Furthermore, it was also found that the subsequent mixing was 
a unique function of (x/x^) . Now even though the and 
injection velocities and densities may be modified, neither 
their mass flow rates nor the density and velocity of the air 
will be affected. Hence, the potential core length is un- 
affected, and the mixing process in the duct will not be 
significantly changed by varying the injection velocities. 


B. DISCRETE TUBE INJECTION 

Practical engineering design requires that a discrete tube 
injection pattern be utilized, rather than the radial (annular) 
configuration studied in the preceding analysis. Calculations 
have been performed which indicate that injection from discrete 
tubes, after a short mixing distance, is equivalent to annular 
injection. These results show that sets of 10 oxidizer and 
fuel injectors located as close as possible to a radial position 
of 2.25 inches (from the center) gives the best (10 percent 

variation about mean temp.) profile. 

The analysis consists of scaling a single fuel (H 2 ) and 
oxidizer (0^) central injector where the propellant flow area 
is enclosed by a radius r .056 ft. Case F shows that 

the edge of the mixing region at L = .8 ft is - *425 ft. 

If now we assume 10 "tubes" having the equivalent fuel and ^ 

oxidizer areas, we get: 


eauiv 


^H2 ^ ^2 
10 


(6.52 X 10 ’^) (37 .5 X 10 S 

10 


= 4.4 X 10 


-4 


ft. 
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The equivalent radius is: 


equiv 


'4.4x10 


I ■■ ■ _2 -2 

1.4x10 = 1.2x10 


equiv 


= .012 ft. 


By scaling# we have: 


^equiv 


where the edge of the mixing region is at 


= 425 =» .091 ft. 

W ^ (.056) 


= 1.1 in. 


“' 2 ? 2 > 


/. . 

mix 


equiv 


where the length of the mixing region is, 


^ 3 i-OlU ^ 3^71 ft. 

W (.056 


= 2.06 in. 

As a typical configuration, let's assume the discrete injector 
tubes are located at a radial position, 

from the center. Then the distance between tubes is approximately. 
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tubes 


^ 2(3.14) (2.5) 


= 1.57 in. 


The interaction radius* obviously is at 


_ C ^ 1.57 
interaction 2 2 


= .78 


In this interaction radius, the interact ionlength is 


interaction niix 


(r ) 

irteraction 

(y • ) 

'mix 


= 2.06 


( 1 . 1 ) 


= 1.46 in. 

Because interaction talces place in such a short length, it 
is concluded that a practical representation (i.e. one that can 
be treated analytically) is the annulus model. Further since 
condition F is the worst mixing case, a design utilizing 10 sets 
of injectors is a conservative one, since better mixing will be 

achlBved for all other conditions* 

Similar results may be obtained by gross scaling. For a 
single central fuel, oxidizer injector, we conclude that a 
length of 15 ft is required (for condition F) to achieve 

uniform mixing. 
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Since the practical length of the combustor is 5 ft, we require 

a scale reduction of 15/S = 3. This implies that the central 
. . 2 

injector area (mass flows) is reduced by (3) =9. Therefore, by 
gross scaling to achieve the total mass flow for the complete 
system, we require 9 sets of injector tubes. 

The analysis for the annular configuration is used to find the 
optimum radial location. Temperature profiles for case F at 
L=5 ft are first plotted for various radial positions of the futl- 
oxidizer annuli, as illustrated in Figure III-18. Next, an area 
pattern factor, A , defined as the area for each curve in 
Figure III-18, that is measured below its maximum horizontal tangent 
point. The area pattern factor for each injector radial position is 
plotted in Figure III-19. It may be observed that the minimum area 
pattern factor occurs at an injector radial position of about 
r = 2.25 inches. This radial location results in an approximate, 
temperature profile of 10 percent variation about mean temperature. 
It is concluded that fuel and oxidizer Injection locations should be 
as close as practical to this radial dimension. 

C. VITIATED HEATER OPERATION AT GASL 

A compilation has been made oi vitiated heater operating data 
as obtained in a GASL combustor which comes closest to the LRC 
design. 

Figure III-20 shows the general arrangement of GASL Combustor 
No. 2. The injection configuration and the ignition sources 

are quite similar to the one designed for the LRC heater . The 
relative locations of the injectors and the ignitor as installed 
in the heater are shown in Figure III-20. Oxygen is introduced 
from an annulus surrounding the igniter; hydrogen injection is from 
twelve tubes located outside the oxygen annulus. The hydrogen in- 
jection direction from the twelve tubes is approximately 45° to the 
combustor horizontal centerline. 
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Operoting data for Combustor No. 2 have been obtained at five 
conditions, •ihose data poihts. which correspond to various test 
programs undertaken with this facility, are indicated in 
Figure III-21. The relationship of these operating conditions to 
the operating envelope and the NASA design points is also provided 
by Figure III-21. Table III-2 gives the pertinent operating data 
corresponding to the five operating conditions. 

Temperature profile data, taken at the test section, with 
cooling air injection at the throat, are presented in Figures III-22 
and III— 23. The profile data correspond to operating conditions 2 
and 3, respectively, in Table III-2. These profiles were adequate 
for our purposes since the test article dimensions are obviously 
always smaller than the test section size, thus permitting locating 
the model hardware in the core region of the flow field where the 
temperature is essentially uniform. Profile data for the other 
operating conditions aro not available* 

In presenting this additional data, our purpose is to show that 
a burner design, operating at conditions within the LRC design 
envelope and similar in design concept to the LRC heater, has been 
suc'.essfully operated. It should be emphasized that point by point 
comparisons are not possible because of differences in design . 
details and in operating conditions. 
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SECTION IV 


ENGINEERING & DESIGN CALCULATIONS 
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A. HEAT TRANSFER 

1. Hgat Transfer to Combustor Wall 

. The combustor walls are subjected to flow conditions .llustratec 
in the performance envelope of Figure IV-1. defined by points k Z D' Z' 


E F G. 

Since efficient mixing and combustion is based on a parallel 
flow configuration, the determination of peak steady state convective 
heat flux can be closely approximated by assuming fully developed tur- 
bulent flow in a tube. Considering the effect of partial dissociation 
of the combustion mixture in the higher temperature range (above 
4500°R) . a heat flux correlation expressed in terms of enthalpy differ 
ence, derived from Reference IV-1, determines the convective heat 

transfer to the coinbustor wall- 

The convective heat flux equation used is derived as follows: 
For large temperature differences between gas and 
wall, the convective heat flux correlation of 
Reference IV-2 is: 


h D 


= .023 


p V D 
f s e 


1.8 


(Pr) 


1/3 


where 

- 2 o„ 

n = film coefficient, BTU/sec-ft - R 
D = hydraulic diameter, ft. 
k = thermal conductivity BTU/sec-ft- R 
p = density, Ib/ft^ 

V = velocity, ft/sec 
M = viscosity, Ib/ft-sec 
Pr = Prandtl No. 

Subscripts 

f = evaluation at film temperature (average between stream 
and wall) 

s » evaluation at stream temperature 
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Since k 


where 


Since 


where 


and 


we may 


where 


Since 




= — ^ — , we can substitute for k, to obtain 
f Pr f 


h = 


.023 (zV )'®(m_)*^(C 
r s t 

, ,2/3 

(D ) (Pr) 
e 


, o 

C = specific heat, BTU/lb- R. 
P 


o 

'C 


h(T -T ) 
s w 


d = convective heat flux, BTU/sec-ft 
^c 

o^ 

T = temperature, R 


Subscript 

w = evaluation at wall 


ff''s 


c V 
s s 




T 



again substitute to obtain 


^c = 


8 T 8 

.023 (i)’ t^)' (Of) -"(Op) 

(D )'^ (Pr)^^^ 
e 


0) = total gas flow rate, Ib/sec 

2 

A = flow area, ft 


AH - enthalpy difference between hot gas and cold wall, BTU/lb 
H = enthalpy of gas at stream temperature conditions, BTU/lb 
H = enthalpy of gas at wall temperature conditions, BTU/lb 
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Thus: 


.023 (f)-® 

A 


^ ^ ) 


.8 


(Mf) 


.2 


(H -H ) 
g w 


(D ) 
e 


.2 


(Pr) 


.67 


In Reference IV-3 and 4, it is recommended that properties 
be evaluated at a reference enthalpy when there is a wide variation 
of gas properties with temperature. This is especially true at 
specification condition Z' where, at 5460 R, the gas has a very high 
specific heat and some dissociation. 

In References IV-3 and 4 the reference enthalpy is evaluated 


as 


= H ^ + .5 (H -H + .22 (H -H . ) 
* st w St r St 


where subscripts 

* = evaluation of reference enthalpy 
st = evaluation at static stream conditions 
w = evaluation at wall conditions 
r = evaluation at a recovery enthalpy condition 


Since the gas velocity is much less than Mach 1, 

H = H ^ (=H ) 
r st g 

and we have 

.22(H -H » 0 

r st 


so that 




.5 (H ) 
st w 

.5 (H +H ) 

9 w 
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In the report we have thus used: 


.023 (i/A)'® (T /T^) *^(H -H ) 

g ♦ * gw 

/r. » -2 X .67 
(D ) (Pr) 
e 


In Reference IV-1 the correlation: 


q = 


.0296 (,*i;/A)*® (m )*^ (H -H ) 

s gw 

iT\ \ *2 /_ \ *67 
(D ) (Pr) 
e 


may be derived from relations presented in that reference, where 
properties are evaluated at stream conditions. Similar results 
would be obtained at condition Z since: 


T .8 

.023 ( ) (u*)* ^ .0296 (u ) 


.2 


An "Enthalpy-Temperature Curve," Figure IV-1, is used to 
determine the difference in enthalpy (:.H) between the gas at stag- 
nation conditions (H ), and the cooled wall (H ). The values of 

g w 

enthalpy are obtained from the output of the same computer program 
used to determine the equilibrium composition of the gaseous system. 

In particular, for specification condition Z' note from the 


curve for 


= .41, that: 


H aj +100 BTU/lb at gas stream temperature T = 5460 R (=5000 F) 

g 9 

and 

H a; -2300 BTU/lb at wall temperature T = 800 R (=340 F) 
w w 
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These absolute values are considered approximate since they are 
plotted for 600 psia initial pressure and 41 percent by weight 
water vapor constant, whereas condition Z' is at 750 psia with a 
slightly lower water vapor content. However, the difference in 
enthalpy (AH) is considered sufficiently accurate to be used for 
heat transfer calculation purposes. 

Convection ; 

■ 023U/A)-°(T^/T.)-°(u.)-^ 

' (D )‘^ (Pr)'^’ 


Sample Calculation for Spe c. Condition 

i = 11.0 Ib/sec 

T = 5000 F = 5460 R 
o 

p 750 psia combustor pressure (ref.) 

H 2: +100 BTU/lb @ T = 5460 R (Equilibrium Data) 

g y ...» 

^ -2300 BTU/lb @ T = 800 R (Equilibrium Data) 

D = 9 in. dia. = .75 ft combustor diameter 

A = .785 (D = .785(.75)^ = .432 ft. 
e 2 

(i/A) = 11.0/. 432 = 25.4 Ib/sec-ft 


* 

H = 


H +H 


w _ 100 - 230^ _ _ BTU/lb 


T* = 3600 R » ar 4.5 X 10 Ib/ft-sec 

(T/T*)= (5460/3600) = 1.52 Pr 2 : .75 
H -H = + 100 - (-2300) * 2400 BTU/lb 

’ ! ,n23(25.4)-^ (] (24001 

q = TTi .67 

(.75) (.75) 


= 179 BTU/sec-ft' 
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Radiation Heat Transfer 

The radiation heat transfer calculations employ values of gas 
emissivity found by using the "upper limit" data of Ferriso, et al. 

(Reference IV-5) , which are determined by experimentation plus 
analysis and extrapolation. 

The heat transfer by radiation can be calculated from the 
simplified equation of Hottel and Egbert found in Reference IV-6: ^ 


q = E 6 [E T 

r o g g 


- a . T ^3 
gb w ^ 


where : 


E = 


E = 


T = 


T = 
w 

'gb " 
6 = 


radiation heat flux from gas to wall, BTU/sec-ft 
effective emissivity coefficient for the wall surface 
emissivity of the gas at temperature T^ 
temperature of the gas, °R 
temperature of the wall, °R 
absorptivity of the gas 

“12 , 2 o 4 

Stephan-Boltzmann constant, .48 x 10 BTU/sec-ft - R 


Because T » T , the term a , T may be considered negligible and we 
gw gb w 

therefore have: 

q - E 6E T ^ . 
r egg 

The values of gas emissivity, E^, were determined by a linear extra- 
polation of the semi-log plot of Reference IV-5, and for the 
condition of interest in this application results in values of the 

order of E = .45 to ,50 for temperature between 4000 to 5000°P, 

9 

For the effective wall surface emissivity coefficient, a value 

of E = .5 is used. This value is considered as a conservative 
e 

estimate between Hottel* s approximation given in Reference IV-6 of: 
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E =-a^ 


E +1 


where E = emissivity of the wall surface considered as a grey 
radiator, (this equation was found to be a fair approximation for 
values of E = .80 to .95, the range most frequently encountered 
in industrial practice), and the other extreme where 


for very low values of emissivity, such as = .072 for ’’comuT.ercial- 

scrapped shiny, but not mirror-like" presented by Hottel in 
Reference IV-7. Since the maximum copper wall temperature wilx not 
exceed 560°F, a value of emissivity of = .57 for a copper plate 

heated at 1110°F is considered too high for this application, and 
therefore, the value of E^ = .5 is considered a good practical 
effective emissivity for design purposes. Radiation heat transfer 
results are plotted for varying combustion temperature at total 
pressures of 750 and 375 psia in Figure IV-2. 

Sample Calculation ; 

Total Pressure, = 750 psia 

T = 3130°K = 5632°R 
T 

Mole Percent of H2O = 50% 

Partial Pressure Water, psia = 25.5 atm 

Characterizing Length, L = .9D = .9(9 in.) =8.1 in. - 20.6 cm. 

PL = (25.5) (20.6) = 525 cm-atm 

From Figure 5 of Reference IV-5, the extrapolated value of = .48 
= (5636)^ = 1000 X 10^^°R^ 

g 

^ • 5 12 2 

/ = E 6E T = (.5) (.48x10"^^) (.48) (1000x10^^) = 115 BTU/sec-ft 

^r egg 
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TOTAL HEAT FLUX 

At specification condition Z', we note that the convective heat 
flux, q = 179 added to a radiation rate of = 96 found from 
Figure IV-2, results in a total heat fl’-ix of = 2 75 BTU/sec-ft . 
Figure IV-3 is a curve of the heat flux to the combustor walls at 
all specification envelope conditions. 

2 . Water Cooling Calculations 

The water cooling calculations are based on the forced con- 
vection heat transfer correlation of Seider and Tate for tubes, 
found in Reference IV-8. 

Calculations to determine wall temperatures as a function of 
heat flux with water velocity as a parameter were perfcrm.ed and 
the results presented in Figures IV-4 and IV-5 for wall gas 
side and water side temperature respectively. 
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Wall Temperature based on Available NASA Water Supply 
The available NASA-LRC water supply of 350 gpm at 480 psi 
corresponds to a water coolant velocity of V = 14 ft/sec in tne 
combustor liner passages. This results in increased uJl tempera- 
tures and thermal growth for any given heat flux. Figure IV-6 
shows the wall temperatures with V = 14 for varying heat flux. 

Also shown are values for V = 20 ft/sec. 

Liner Temperature Response at Various Thermocouple Depths 
An analysis has been made to determine the zirconium copper 
liner wall transient temperature for various fractional wall 
depths as a parameter. The calculations were performed assuming 
the wall to be a plate with sudden exposure to a constant heat in- 
put on one side with the other side insulated (no water cooling) , 
to simulate a malfunction type condition. The temperature response 
data of Reference IV-12 was used, where: 





where : 
F 

o 

T 

a 

8 

6 

k 

• 

q 

t 

t 


p 

c 


= Fourier No. = ^ 

6 

- dimensionless temperature , 

k ft^ 

= thermal diffusivity of body material = — « — ~ 

pc sec 

" time, sec 
= plate thickness, ft. 

= thermal conductivity of body material, BTU/sec-ft-°F 

2 

= rate of heat input, BTU/sec-ft 
» temperature, °P 

o ' 

B initial temperature, F 

. 3 

<■ depth of body material, Ib/ft 

specific heat of body material, BTU/lb-^F 


P 


X - T 


X 


depth ratio 
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X = depth, ft. 
(x = 0 at X = 0) 
(x = 6 at X = 1) 


For zirconium copper 

K = 200 BTU/lb-ft-°F 

= 5.56 BTU/sec-ft-F 

p = 559 Ib/ft^ 

C = .0915 BTU/lb-°F 
F 

a = = 1.09 X 10 ^ ft^/sec 

oC 

P 

6 = .5/12 = .417 X lO"^ ft 
SAMPLE CALCUIATION ; 

2 

At 9 = 1 sec, q = 200 BTU/dec-ft 

p _ ci _ (1.09x10-3) (1) 

'o ■ 2 “ (.174 X 10“^) 

P 

= .627 

At X = 0, X = 0 

From Chart 43 of Reference IV- 12 

T = .96 
t-t» = * T q 

(.417xl0~^) (.96) (200) 

_ ' 10-i^) ^ 

= 144® AF 

t = 80P 
o 

t = 80 + 144 = 224 °P at e ■ 1 sec 
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The results of these computations are shown in Figures IV-7 and 8. 
Figure IV-7 shows the temperature response at one second as a 
function of heat flux, with depth as parameter, and Figure IV-8 
shows the temperature response at q = 275 as a function of time, 
again with depth as parameter - 

From Figure IV-7 vre see that at a depth of X = 2/3, x = .333 
in (distance of .167 inches from insulated side) there is a 
temperature rise of 110 after 1 second exposure. This 

compares to a rise of 200 ^F at X = 0 (hot surface) . A thermo- 
coup.’.e embedded a distance of .167 inches from the cold side 
(X = 2/3) is therefore considered adequate. 

THl-:R^tOCOUPLE INSTALLATION 

To secure the 1/16” diameter sheathed thermocouple in place, 
a thermally conductive cement is to be used. This is preferred 
rather than a true soldering operation in order to avoid heating 
the entire liner after machining. 

The cement specified is "Eccobond Solder 57C,” an epoxide 
based room temperature curing conductive adhesive. It is used in 
applications where conventional hot soldering is impractical. 

The thermal conductivity given by the manufacturer (Emerson and 
Cuming, Inc.) is k > 200 BTU/ft^/hr/°F/in. 

3 . Air Cooling Analysis 

An analysis was performed to determine the suitability 
of using heater supply air as a combustor wall coolant, prior 
to injection into the main combustion chamber. 

The analysis, evaluated at condition Z*, indicated that 


air was not an appropriate coolant choice for this application. . 

The analysis follows. I 

!■ \ 
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Radial Expansion of the Combustor Liner 
Analyses have been made to determine the radial thermal 
expansion of the liner and the corresponding clearances required. 
Radial expansion results from heating of the copper liner from 
80 F initial temperature to a final average wall temperature. 

At specification condition Z', with nominal water coolant 

2 

passage velocity, V = 14 ft/sec, and heat flux, q = 275 BTU/sec-ft , 

we have an approximate average copper temperature rise of aT=365 :.F. 

Radial Expansion of Liner at Instrumentation 
Port Locations 

Access ports for instrumentation are located at Section B-B 
as shown on GASL drawing HE 1702 (NASA drawing No. L522238) . A 
clearance space must be provided between the inner and outer liners 
to accomodate thermal radial expansion of the inner liner. The 
average temperature rise of AT == 365 aF was used to determine the 
expansion. 

The diametral radial thermal expansion is approximated by; 
where ; 

6, = diametral growth, inches 

d Q 

a = coefficient of thermal expansion, in/in- F 

o 

AT, = average bulk tem.perature rise of copper liner, F 
b 

D s diameter of liner, inches 
» (9.3 X lo”^) (365) (10.0) 

a 

= .034 inches (diametral growth) 

The design will provide for .040 inch diametral clearance at 
Section B-B. 

Radial E.iPan:jion of Liner Downstream Flange 
The liner downstream sealing flange is shown on GASL draw- 
ing HE 1700 (NASA drawing No. L-522236) . 

To determine the radial expansion at the downstream sealing 
flange, the temperature profile through the flange must be known. 
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A graphical method found in Reference IV-10, was used to 

approximate (for design purposes) the temperature profile. The 

method requires sketching into the flange regxon a net of 

isothermal-adiabatic lines, carrying out the construction by 

trial and error until the criterion that these lines meet at 

right angles at all points of intersection is satisfied. 

The results of this approximation are shown in Figure IV-10 

for a heat flux to the liner wall assumed constant at 

2 

q = 275 BTU/sec-ft and with the water flow area restricted so 
as to obtain a coolant velocity of = 40 ft/sec in this region. 
The coolant film coefficient in this curved-entrance type flow 
region is approximated from the method of Jeschke for helical 
flow found in Reference IV-11: 


iLl 

h 


1 + 3.5 



where: 

h'/h = ratio of coolant heat transfer coefficient in 

helical flow (curved) to that in a straight pipe 

D = pipe diameter (equivalent diameter, D^, used) , ir>ches 

D.. = helix diameter (2 x corner radius used) , inches 

He 


To obtain V = 40 ft /sec, coolant passage must be reduced to: 


t = .25 (^) = .088 in 

D = 2 (.088)* .176 in 

e 


D„ * 2 X .25 
Ha 


.50 


h' 1 L n • M76) , , 1 

— * 1 + 3,5 ;■ ■ 'v ** 1 + 1.23 
h (.50) 


2.0 
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For the value h'/h c_. 2, we can use the water coding character- 
istic curves of Figure IV-4 for V = 80 to obtain approximate 


wall temperatures at q = 275, V 80 


T 

WALL 

COLO SIDE 


= 160®F 


T 

WALL 
HOT SIDE 

T 

WALL 

AVERAGE 


= 360°F 


l60-»-360 

2 


= 260°F 


Observation of Figure IV-10 indicates that an isotherm of 
o 

T = 260 F passes through the flange (and O-ring) section, and 
this value is used as the average temperature for determining 
a copper bulk temperature rise. 

= 260 - 80 = 180 L^F 


the radial expansion is approximated as 

C-d = a ATgD 

= (9.3 X 10“^) (180) (11.0) 
s .018 inches (diametral) 


Therefore, in order to provide for radial expansion at the 
liner downstream sealing flange, a radial clearance (diametral) 
of 0.018 can be provided. Static 0-ring sealing can be attained 
with this clearance. 
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The pressure drop at the entry flange, where V = 40, is due 
to a contraction, = *5, two right angles, = 1.0, = 1.0, 

and an expansion = 1.0, so that 

AP„, = % 

Flange • 2g 


where K = pressure loss factor 


AP, 


Flange 


= [.5 + 1.0 + 1.0 + l.Oj [ 


(62.4) (40) 

9.28 X 10^ 


■] 


= (3.5) (10.8) 

= 40 psi 

The pressure drop in the boss region, with a slot height of 
.188, consists of a contraction loss, a friction loss and an 
expansion loss. For a water capacity of 350 gpm, the water 
velocity past the bosses is 
w 


V = 


pA 


w 


where : 


A = water flow area between bosses (ft ) 
w 

A - (irD) blocicage length) t 


w 


= f (3.14) (10.188) - [(4)(1) + (2) (2.5)]] (.188) 

= (32.0-9) (.188) 

*>4.33 in^ = .030 ft2 


Therefore: 
V < 


/ 152, ) ( — L_ ) ^ — L. 

' 7.2 ^ ' 62.4 * ' .030 
26.9 ft/sec 


) 


107 


tacm <»n.HD fctmcc ftow*Toi>m me. 


•— 7%wbti<i«fT 0t 




OfQuordt 


The boss region entry plus exit losses are 




pV^ 


(Ks+Ke) 



(.5+1.0) 


, (62.4) (26.0)^ 

9.28x10^ 


= (1.5) (4.55) = 6.83 psi 


The frictional loss in the boss area is: 

■ p , f, i , ££ 

D ^ 2g 


= (.03) ( 


.188 
= .48 (4.55) 

= 2.2 psi 


) (4.55) 


The pressure drop through the boss area is then 




The pressure drop at the liner exit consists of a right angle 
= 1 and an enlargement = 1 at V = 14 ft/sec 




''“exit = ^ ■' 


) 


(1+1) (1.32) 


« 2.6 psi 

The total liner pressure drop is therefore: 

^^total * ^^entry ^ ^^length ^ ^^boss ^ ^^exit 
liner flange friction region 

» 40.0 + 4.3 + 9.0 + 2.6 

- 55.9 » 56 psi 
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